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Abstract 

. We review the literature concerning the Hardy inequality for regions in 

Euclidean space and in manifolds, concentrating on the best constants. We 
also give applications of these inequalities to boundary decay and spectral 
approximation. 
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1 Introduction 



Let if be a non- negative second order elliptic operator acting in L 2 (U) subject 
to Dirichlet boundary conditions, where U is a region in or in a Riemannian 
manifold. Also let d be a positive function on U which is continuous and satisfies 
|Vrf| < 1. Traditionally one takes d(x) to be the distance of x £ U from the 
boundary dU, but another possibility is that d(x) is the distance from any closed 
subset of M\U if U is embedded in some larger Riemannian manifold M. 

We say that H satisfies a weak Hardy inequality with respect to d if there exists a 
constant c > and a constant a > such that 

i/i 



|ML< c2(w) + 



is valid for all / G C%°(U), and hence for all / in the domain of the quadratic 
form Q of H. The infimum of all possible c in (IJ) is then called the weak Hardy 
constant. We say that H satisfies a strong Hardy inequality if ([!]) holds with a = 0, 
in which case the minimum possible c is called the strong Hardy constant. 

There are also LP and higher order analogues of the above notion, which we mention 
briefly later in this review. 

In section 2 we describe the method of geodesic integrals for proving Hardy in- 
equalities in higher dimensions. Section 3 describes a method ultimately due to 
Jacobi, while Section 4 gives various miscellaneous results. 
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We then turn to the applications of the HI to the proof of boundary decay. It 
was shown in |2l] that Hardy's inequality can be used to prove the L 2 boundary 
decay of eigenfunctions without any further assumptions. This in turn leads to 
the possibility of controlling the rate of convergence of the eigenvalues when the 
region U is approximated by a family U £ of slightly smaller regions. Very recently 
progress has been made on this problem, |24[] , and we are able to announce bounds 



on the rate of convergence which are sharp in a certain sense. 

Our main results on boundary decay, Theorems 11 and 12, may be regarded as L 2 
analogues of much stronger pointwise bounds on eigenfunctions given in |13, |34 



Note however that our bounds depend only on the validity of (|1|), hold for all 
functions in the domains of the operators, not just for eigenfunctions, and have 
rather precise constants. 

If we abandon interest in the precise value of the constant, and choose d to be 
the Euclidean distance from an arbitrary point of U, then it may be seen that our 
results are related to Morrey space estimates. These have been of considerable 
importance in the theory of elliptic operators, and recently in the proof of heat 
kernel bounds, and we refer the reader to P, |52| for further details. 



2 Geodesic integrals 

The first method which we describe depends upon the one-dimensional case, which 
is the only one Hardy actually studied. We refer to ]I0j for an exhaustive study, 
which involves generalizations to the variable coefficient case of the original formula 

r j/MJ! dx < 4 f°°\f {x) \2 dx 

Jo x A Jo 

valid for all / G C c °°(0, oo), and hence for all / G W 1,2 {0, oo). 

Let H := -A^/r in the Hilbert space L 2 {U) where U is a bounded region in R^. 
For every unit vector u G S N ~ l and x G U let 

d u (x) := min{|t| : x + tu ^ U} 

if the set of such t is non-empty, and put d u (x) := +oo otherwise. We define the 
(harmonic) mean distance of x from dU by 

m ( x )- 2 : = l^" 1 !" 1 J i d u (x)~ 2 dS(u). (2) 

It is easy to prove that d(x) < m(x) for all x G U. 
Lemma 1 We have 

N <H 



Am 2 
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in the sense of quadratic forms. If Ai is the smallest eigenvalue of H then 

N 

where the quasi-inradius n of U is defined by 

(j, := sup{m(x) : x G U}. 



Proof See @ or 0, Th. 1.5.3]. 

Applications of the above lemma depend on making assumptions on U which enable 
one to bound m(x) above by some multiple of d(x). The first of these is folklore 



and seems not to have been written down explicitly until very recently; see |p6| , |37 
and the next section for alternative proofs. 



Theorem 2 If U is a convex subset of H N then 



^ 4H 



in the sense of quadratic forms. 



Proof If a is the point of dll closest to x then we can obtain the relevant upper 
bound of m(x) by computing an appropriate integral over the supporting hyper- 
plane at a. See p3], Exercise 5.7]. 



The following lemma is typical of a variety of methods of obtaining crude upper 
bounds on m(x). The hypothesis is valid not only for regions with Lipschitz bound- 
aries, but also for a variety of regions with fractal boundaries, such as the Koch 
snowflake region in R 2 . 

Lemma 3 Suppose that there is a constant k such that for each a G dll and each 
a > there exists a ball B disjoint from U with centre b and radius (3 > ka, where 
\b — a\ = a. Then there exists constants cq, c± such that m{x) < cod(x) and hence 

in the sense of quadratic forms. 



Proof See |TS|, [0, Th. 1.5.4] and |T} Th. 3]. 

The condition of Lemma 3 is not satisfied for regions satisfying a uniform exterior 
power-like cusp condition. In such cases one may prove a modified Hardy inequality 
using Lemma 1, namely 

/^<c 2 (w)+«imi 2 ) (3) 

for some < 7 < 2; see [|28|, p 369]. See also p6|, Th. 3.2, 3.3] where a similar 
situation arises for locally Euclidean manifolds with fractal boundaries. 
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A procedure closely related to the idea of this section was developed for regions in 



Riemannian manifolds independently by Croke and Derdzinski, H], and Donnelly. 



30| . The integrals over straight lines were replaced by integrals over geodesies, so 
the formulation involves the geodesic flow on the unit sphere bundle of the manifold. 
However, both papers are concerned with obtaining lower bounds on the bottom 
eigenvalue, much as in Lemma 1, rather than Hardy's inequality. 

We mention in passing that there is no requirement that one should assign equal 
weights to every direction in Euclidean space. In some cases one obtains a better 
constant in the Hardy inequality by taking an average over a few directions which 
are well adapted to the region in question. 



3 The Classical Method 

The following method goes back to Jacobi, and was used by Barta and Kasue to 
obtain lower bounds on the first eigenvalue, [f|], |33|. It is the easy half of a theorem 



of Allegretto, Moss and Piepenbrink characterising the bottom of the spectrum of 
a Schrodinger operator in terms of the existence of positive distributional solutions 



of the eigenvalue equation, [15, p. 23]. Assume that 

df 




Hf{x) := -2^^T W 



dXn 



where a(x) is a non-negative C 1 real symmetric matrix- valued function and / G 
C 2 (U). Then if is a non-negative symmetric operator and we can use the same 
symbol to denote its Friedrichs extension. 

Lemma 4 Let <fi be a positive C 2 function on U and let V be a continuous function 
on U such that 



Then we have 

in the sense of quadratic forms. 



H>V 



Proof See 0, Th. 4.2.1]. 



The conditions of the above lemma can be weakened to allow a distributional 
inequality. 

Second proof of Theorem 2 If we put <fi := d 1 ^ 2 and use the fact that Ad < for 
any convex set U then the result follows immediately from the last lemma. 

The method of this section can be extended to Riemannian manifolds without 



difficulty We refer to |T0, E5L |1J for a variety of Hardy and Rellich type inequalities 



with explicit constants in Riemannian manifolds obtained in this manner. The 
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following theorem is only one of a range of related results due to Brezis and Marcus, 
|| . In particular they find explicit bounds on the minimum possible negative value 
of a in the theorem when U is convex. 

Theorem 5 If U C is bounded with a C 2 boundary and H := —Adir in 
L 2 (U) then there exists a G R such that 

d- 2 < A{H + a) (4) 

in the sense of quadratic forms. If U is convex then (0j holds for certain a < 0. 

Proof Let <p be a positive C 2 function on U such that (p(x) = d(x) 1 ^ 2 — d(x) for 
all x close enough to dll. The first statement of the theorem follows by applying 
Lemma 4 to </>. 

There are various other improvements of the strong Hardy inequality of which we 
mention just two. For a definitive treatment of the one-dimensional theory see ED 



Theorem 6 If U := {x G R : x^ > 0} where N > 1 then 



u | x 2 N + 4:X N (x 2 N + 



for all f<=C™{U). 

Proof See |8|, Sect. 2.1.6]. 

Theorem 7 If U := (0, a) i/ien 

« 2 m 2 



/ 2 , , 2 dx<4/ |f| 2 dx 
Ju x z (a — xr Ju 



for all feC?(U). 

Proof Put <6 fx) := x 1//2 fa — x) 1 ^ 2 in Lemma 4. 



4 Capacity-based methods 

In this section we mention a few of the very general theorems which involve the use 
of capacity arguments. These have been developed in an L p context, but we only 
treat the case p = 2. If K is a compact subset of U C H N we define its relative 
capacity by 

capfiT, U) := inf | V/| 2 : / G C7f>(tf) and f\ K > l} . 

It is particularly appropriate in this conference to mention one version of the most 
quantitatively precise theorems of this type, due to Professor Maz'ya. 
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Theorem 8 If // is a positive measure on U and 

for all compact subsets K of U , then 

[ |/| 2 d/i < 4(3 f |V/| 2 
Ju Ju 

for all f G C%°(U). Conversely the second inequality implies 

for all compact subsets K ofU. 
Proof See @, p. 113]. 

Our next results are taken from a paper of Ancona, [l[] . We say that U is uniformly 
A-regular if for all x G dU and all r > the harmonic measure w of U D dB(x, r) 
in U fl B(x, r) satisfies w <1 — j3 on U H dB(x,r/2), for some constant (3 G (0, 1) 
independent of x, r. If iV > 3 this is equivalent to the uniform capacitary density 
condition that there exists a constant a > such that 

c&-p(B(x,r)\U) > ar N ~ 2 

for all x G dU and all r > 0. 

Theorem 9 If N > 2 and U C H N is uniformly A-regular then U satisfies a 
strong Hardy inequality with respect to the Laplace operator. If N = 2 then the 
converse is also true. 

Although UH does not provide sharp information about the size of the strong Hardy 
constant, it contains many more results than we have indicated above. An LP 



converse of Theorem |9] for N = p > 2 may be found in |35| , using an appropriate 
L p Riesz capacity. 



5 Miscellaneous results 



The weak Hardy constant c as defined in Section 1 was proved in |22] to be local in 
the sense that it is the maximum value of a certain upper semi-continuous function 
on the boundary, whose value at each point depends only on the geometry of 
the boundary around that point. Various methods of evaluating this function at 
different types of boundary point are described in . 

For the remainder of this section we assume that H := — A^jr. The strong Hardy 
constant is a global invariant of U. It equals 2 for any convex set, but the condition 
of convexity is not necessary for this conclusion. Let 

Up := {re 10 : < r < 1 and < 9 < /?}. 
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Then Up has strong Hardy constant 2 if and only if the internal angle (3 is less 
than or equal to a certain critical value (3 C ~ 4.856 radians, [22|. For larger /3 the 
strong and weak Hardy constants are larger than 2. Similar conclusions hold for 
other plane regions with piecewise smooth boundaries. 

If U is a simply connected region in R 2 then U has strong Hardy constant at most 4 
by 0, |17], Th. 1.5.10]. The proof of this result depends upon a fact from analytic 



function theory, namely Koebe's one-quarter theorem. 

There is an interesting connection between the possible constants in the strong 
Hardy inequality and the Minkowski dimension of the boundary, [^7|. In two 



dimensions there is also a relationship with hyperbolic geometry, which we do not 
pursue. We say that the boundary dU has interior Minkowski dimension a if there 
exist positive constants k\ and k<i such that 



ki£ N~a <\{ xeU: dist{x,dU) <e}\< k 2 e 



N-a 



for all e > 0. The following theorem is adapted from |27], Th. 3.3]. We allow 
a < N — 1 because the theorem is applicable in manifolds, for example if U is 
obtained by removing a compact set K from a sphere endowed with the standard 
metric. 



Theorem 10 If dU has interior Minkowski dimension a > N - 
Hardy constant of U with respect to the Laplacian satisfies 

c(2 + a - N) > 2. 



2, then the strong 



In most of the above lemmas we have restricted attention to Hardy inequalities in 
L 2 . In fact many of the results have been extended to LP with sharp constants; see 



[36], [37| for the proofs of the following two theorems. 



Theorem 11 Let 



inf 



Su\f/d\> 



: / G Wp(U) 



where 1 < p < oo. If dU is smooth then c > p/(p — 1). If in addition p = 2 then 
c > 2 if and only if the infimum is achieved by some f G Wq' 2 (U). 



Theorem 12 If U is a convex set in R and 1 < p < oo then 



u dP 



< 



V 



p — 1 



|v/p 



for all f G W^{U) 



We refer to 0, g |j| H, ||, 0, H || for further L p results, since they do not 
yet have such direct consequences for spectral theory. We refer to [pB], |4~l~f for the 
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analogues for higher order operators, known as Rellich inequalities, and to [|T^] for 
analogues in Orlicz spaces. 



We describe some trace inequalities in JT9| which may be proved using Hardy's 
inequality. We assume that H := — &.dir acting in L 2 (U) where U is a region in 
R w . The theorems are only of interest when U has infinite volume. 

Theorem 13 We have 

tr[e- m ] < {2id)~ N ' 2 f e- m W*) 2 d * x 

for all t > 0, where m is defined by (fjj. 

Theorem 14 If U satisfies the regularity condition 

d(x) < m(x) < bd(x) 

for all x G U then 

2- N (2irt)- N / 2 ( e - 8 ^ NH / d ^ 2 d N x < tr[e- Ht ] < (2id)~ N ' 2 I e~ Nt ^^d N x 
Ju Ju 



for all t > 0. Hence 

for all t > if and only if 

for all t > 0. 



tr[e" m l < oo 



e -t/d(x)* d N x < ^ 
U 



6 Boundary estimates 

The size of the constant s in an inequality of the form 

\f\ 2 



u d s 



< oo (5) 



conveys information about the behaviour of the function / near the boundary of U. 
We conjecture that it is not possible to have s > 2 in the inequality for any region U 
if we only assume that / G Dom(Q), where Q is the quadratic form associated with 
a uniformly elliptic second order operator H acting in L 2 (U) subject to Dirichlet 
boundary conditions. However, if we make stronger assumptions on / then one 
may be able to prove (||) for a larger value of s. The first paper with results 
of this type was PH , where it was assumed that / was an eigenfunction of H. 
Subsequently |yj obtained better bounds for all / G Dom(if), assuming only the 



Hardy inequality. 



8 



Although we have concentrated on L 2 boundary estimates, there is a substantial 
literature on pointwise decay of eigenfunctions and their gradients at the boundary 
Bounds of the type 

\4>n(x)\ < C n <pi(x) (6) 



are immediate consequences of intrinsic ultracontractivity (IU), |2l| 17], in which 
a major ingredient of the proof is the existence of an inequality 

<j>x{x) > ad(x) a (7) 

for some positive constants a and a. The proof of (|7J) depends in turn upon the 
Harnack inequality and a boundary accessibility property. The BAP was proved 



in [^8|, [TTJ for Lipschitz domains, but Ancona and Simon commented that it holds 
under a suitable twisted interior cone condition, i.e. for John domains, [^, p 98]. 
Finally Banuelos gave a detailed analysis of the relationship between (g), IU, John 
domains, Holder domains, NTA domains, etc. in Q. 

Pointwise bounds on the gradients of the eigenfunctions <fi n of —Adir and of 
Schrodinger operators with potentials in restricted Kato classes acting in L 2 (U) 
are proved in [11], [13|, |34], |], || in steadily increasing generality. The best upper 



bound is for IU domains and is in ||, while the best lower bound is for Lipschitz 
domains and is in |J. The inequalities are of the form 

\V<p n (x)\ < c n (t)i{x) / d(x) 
|V0i(x)| > c X (j) 1 {x)/d{x), 

the latter being for x close enough to the boundary. 
One may also obtain upper bounds of the form 

\(t>n{x)\ < C n d{xf 

for explicit but non-optimal constants c n , /3 which depend only on the eigenvalue 
A n , the dimension and the constant a in the uniform capacitary density inequality, 
|| . For an open simply connected region in R 2 the bound 

\<p n (x)\ <c n d(xy/ 2 



is proved in [|, §, £|; in this case the power 1/2 is sharp. 

We finally present some new results on L 2 boundary decay, taken from |24] . Let U 
be a bounded region in R^ and let H := —Adir acting in L 2 {U). Let d(x) denote 
the distance of x from some closed subset of R \U. We make no assumptions on 
the boundary dll apart from the validity of (^) for certain values of c > 2 and 
a > 0. We are then able to draw the following conclusions about the boundary 
decay of functions in the domain of H. We have proved in that the powers of 
e in these theorems are sharp, and conjecture that the constant cq is also sharp. 
Analogues of the theorems for uniformly elliptic operators in divergence form are 
proved in p4|. 
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Theorem 15 If f G Dom(if) and e > i/ien 

/ l/I^Cde^^lK^ + oJ/HalK^ + a^/lb 

7{s:d(a;)<£} 

w/iere 

,2+2/c 



c := c 

Theorem 16 If f G Dom(iJ) and £ > i/ien 

/ IV/l^d^Kff + o^lhlK^ + a^Vlb- 

J{a;:(i(x)<£} 

:= c 2 / c + c 2 / c (l + c) 2+2 / c . 



Ci : 

Corollary 17 //if/ = A/ /or some A > 0, ||/|| 2 = 1 and e > t/ien 

/ |/| 2 <c e 2+2/c (A + a) 1+1 / c 

J{x:d(x)<e} 

and 



{x:d(x)<e} 



|V/| 2 < Cl e 2 / c (A + a 



,1+l/c 



Corollary 18 If U is a simply connected proper subregion of R 2 then 

I „ |/| 2 <32 £ 5 / 2 ||iJ/|| 2 ||i/ 1 /4 / || 2 



and 



1 



' {x:d(a;)<e} 

/or a// / G Dom(if) and £ > 0. 



|V/| 2 < lUe^\\Hf\\ 2 \\H^fl 



Proof We put c = 4 and a = in Theorems [15| and [16 . 

We use the results above to consider the effect on the spectrum of H := — ^dir of 
replacing the bounded region U by a slightly smaller region U e such that 

{x G rj : > £} C C/ £ C 17. 

If X n (U e ) denote the eigenvalues of the operator H e defined by restricting H to 
L 2 (U £ ) where we again impose Dirichlet boundary conditions, then variational 
arguments imply that X n (U) < X n (U E ) for all n and e > 0. Our theorem below 
provides quantitative estimates of the difference, again only assuming ([[]). The 
first version in |H[] did not obtain what we believe to be the sharp power of e 
given below. Pang, obtained the result of Theorem ^ for n = 1 for simply 
connected plane regions by a method involving conformal mappings, improving his 
own earlier results in ||42|| . See for a more general version of the theorem below, 
and its proof. 
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Theorem 19 There exist constants c n for all positive integers n such that 

X n (U)<X n (U £ )<X n (U)+c n e 2 / c . 

We finally mention that there is extensive literature which compares X n (U) with 
X n (U\K), where K is a compact subset of U which has a small capacity in a suitable 
sense; we believe that J3) is one of the earliest contributions to this subject, often 
known as the crushed ice problem. See [2^] for a survey, including an explicit 



asymptotic formula for the difference of the eigenvalues in the limit of small Cap(fT) 
and also estimates of the difference for n = 1, both proved in the abstract setting of 
regular Dirichlet forms. See also , where estimates of the difference for n = 1 in 



terms of an appropriate definition of capacity are obtained in an abstract contect 
applicable to higher order elliptic operators. 
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